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Characteristic properties of magnetic field efi'ects on spin selective geminate and bulk electron-hole 
polaron pair (PP) recombination are analyzed in detail within the approach based on the stochastic 
Liouville equation. Simple expressions for the magnetic field (B) dependence of recombination yield 
and rate are derived within two models of relative PP motion: free diffusion and diffusion in the 
presence of well (cage). The spin evolution of PPs is described taking in account the relaxation 
induced by hyperfine interaction, anisotropic part of the Zeeman interaction induced, as well as 
Ag-mechanism. A large variety of the B-dependences of the recombination yield Y{B) and rate 
K{B) is obtained depending on the relative weights of above-mentioned mechanisms. The proposed 
general method and derived particular formulas are shown to be quite useful for the analysis of 
recent experimental results. 
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I. INTRODUCTION 

Magnetic field effects (MFEs) on various processes in 
organic semiconductors are actively studied for many 
yearsii"— These studies concern different types of MFEs: 
the magnetic field dependent photoconductivity;^ and 
photoluminescence^i^, magnetoelectroluminescence^^ii^ 
magnetoresistance^"— etc. 

The mechanisms of MFEs in these processes are the 
subjects of hot debates for many years. It is, however, 
widely acceptedfi^"— that a large number of MFEs result 
from the effect of the magnetic field (B) on the spin se- 
lective reactions with participation of paramagnetic par- 
ticles: polarons (P) and triplet excitons (T). In a large 
number of them the key stage is the recombination of 
pairs of electron (e) and hole (h) polarons, i.e. particles 
with electronic spin 1/2^^1^^^^^ which are called here- 
after polaron pairs (PP). The T — P quenching T — T 
annihilation are also believed be the important spin se- 
lective process which can give significant contribution to 
the observed MFEs in organic semiconductors f^S^"— 

There are some reviews of experimental and theoretical 
works on MFEs in organic semiconductorsj^^"— Recent 
extensive experimental investigations, however, inspire 
further theoretical studies of MFEsi^i^ii^iii Despite ev- 
ident progress in these studies there are still many prob- 
lems to be discussed. 

In particular, close attention has been attracted to the 
e—h PP recombination mechanism of MFEs in disordered 
semiconductors!^ ' "'^''i The important problem consists in 
proper treatment of the effect of polaron migration and 
disorder of spin dependent interactions, giving rise to the 
MFEs. The majority of theoretical works are mainly 
based either on numerical or somewhat simplified analyti- 
cal description of the spin/space evolution of PPs^i^iilii^ 
though fairly high accuracy of above mentioned recent 
measurements motivates more detailed theoretical inves- 
tigations, which could allow for obtaining sufficiently ac- 



curate and rigorous formulas for the MFEs. 

In this work we discuss the method of describing 
specific features of the PP recombination mechanism 
of MFEs in disordered organic semiconductors. The 
method is based on the diffusion approximation for hop- 
ping polaron migration, fairly reasonable at long times 
of MFE formation. The hopping motion is assumed to 
lead not only to stochastic spatial evolution but also to 
fluctuations of disordered spin dependent hyperfine inter- 
action (HFI) and anisotropic part of the Zeeman interac- 
tion (AZI), resulting from the anisotropy of 17- factors of 
polarons. The fluctuating HFI and AZI give rise to spin 
relaxation, which is described by Bloch-type equations, 
valid in the realistic limit of hopping rates much larger 
than these interactions (in frequency units). 

The kinetics of MFE generation is determined by the 
spin/space evolution of PPs, which is described by the 
PP spin density matrix. In the above-formulated ap- 
proximations this matrix satisfies the stochastic Liouville 
equation (SLE),i^'^>25 

With the SLE approach we analyze the properties of 
MFE for geminate and bulk processes within two mod- 
els of PP relative motion: free diffusion and diffusion in 
the presence of a potential wel h^^'^'' The PP spin evo- 
lution is described taking into consideration the above- 
mentioned HFI and AZI induced relaxation, as well Ag- 
mechanismi^ In these two models of relative motion sim- 
ple expressions are derived for the MFE, i.e. for B de- 
pendent recombination yield y{B) and rate K{B) (in 
geminate and bulk precesses, respectively). 

The analysis with obtained formulas reveals different 
types of Y{B) and K{B) behavior, as B increased, for 
HFI and AZI induced relaxation mechanisms: decreas- 
ing and decreasing, respectively. The combination of 
these mechanisms is found to result in a large variety of 
non-monotonic Y{B) and K{B) dependences. The extra 
contribution of Ay-mcchanism can lead to some addi- 
tional specific features of MFEs behavior at large mag- 
netic fields B whose specific features appear to depend 
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on the mechanism of relative motion of polarons. 

In our discussion we also concern some possible appli- 
cations of obtained expressions to the interpretation of 
some recent experimental results. 



II. MECHANISM OF POLARON MIGRATION 

The characteristic properties of MFEs on the migra- 
tion assisted PP recombination are, naturally, essentially 
determined by the mechanism of polaron migration. 

There are a number of models of migration in dis- 
order semiconductors. One of the most popular is the 
Miller-Abrahams modeL ^^i^^ In this model the rate Wij 
of hopping from the site j to the site i is written as 
Wij = Wj^,^_ = woe"^^^'"-^^")^-^*"-^^)/^''^^), where wq is 
the characteristic rate constant, 9{x) is the Heaviside step 
function, is the distance between these sites, and Ej 
and Ei are the energies of initial and final states, respec- 
tively, assumed to be randomly distributed parameters 
whose (broad) distribution functions are determined by 
a number of intra and interpolaron interaction.— 

Even in this relatively simple model the kinetics of the 
space/time evolution of charge carriers can be obtained 
only numerically. 

It is worth noting, however, that the MFEs are de- 
termined by very long times of order of the character- 
istic time Tg ~ 10~^s~^ of spin evolution of e- and h- 
polaronS)^ which is much longer than the average hop- 
ping time Wq^: Wq^ (see Sec. VIII). At times 
t > Tg ^ Wq^ in the wide region of parameters of the 
model the space/time evolution is fairly reasonably de- 
scribed by the diffusion approximation. One of indica- 
tions of this fact is a reasonably good accuracy of the 
Langevin formula for the recombination ratei^ 

In the diffusion approximation the migration kinetics 
of z/-polaron (v = e,h) is described by the time dependent 
probability distribution function pi,{r^,t) in the contin- 
uum space {tc^}, obeying the Smoluchowski equation 

In ea. (|2.ip ^ ^I^q is the effective diffusion coefhcient 
for z/-polaron, in which and wq^ are the characteristic 
hopping length and hopping rate. Note that eq. (12. ip 
is actually a selfconsistent equation, which incorporates 
many particle effects, showing themselves in possible de- 
pendence of the diffusion coefficient on the concentration 
of polarons, as well as concentration and temperature 
dependence of the effective potential u{r^)^ 

Our further analysis of MFEs will be based on the dif- 
fusion approximation (|2.ip . Nevertheless, some specific 
features of hopping kinetics at short times, predicted by 
the Miller- Abrahams model, appear to be important as 
well for the description of MFEs (see below). 



III. SPIN EVOLUTION 

Quantum (spin) PP evolution plays the key role in the 
generation of MFEs on PP recombination. PP recombi- 
nation is known to be a spin selective process with the 
rate depending on the total electron spin of the pair of 
e- and /i-polarons: S = Se -I- S/i, where 8,^ is the spin of 
I'-polaron (v — e,h). In our discussion we will assume 
that the rate is non-zero only in the singlet (S) state 
corresponding to the total PP spin S = 0.^*'^^ 

The spin evolution of e- and /i-polarons is described 
by spin density matrices p^, (v = e, h), in the two state 
Hilbert spaces, in which two states \i>±) correspond to 
two spin projections onto z-axis, taken to be parallel 
to the vector B of external magnetic field: S,yji'±) = 
±i|j/-i-). In general, these matrices can be represented in 
terms of expansion in bilinear combinations of the states 

W^^') = Wt.){'^t.'\- (3.1) 

It is important to note that can be considered 

as vectors in some space, called hereafter the Liouville 
space and denoted as {i^^^'}, in which the conjugated 
vectors are defined by the relation {v^_^^'^ \v^^^'^) = 

S^_^^^6f^>^^i^. Noteworthy is also that this definition im- 
plies linear independence of the vectors Ivfi^.') and jt'/i'/i), 
i.e. (v^'^lv^^i) = 0. In {j/^^'}-space, the density matrices 
= J2 , j_Py / \ can be written as a vectors: 

Pu = \Pu) = (3-2) 

In addition to vectors, in the Liouville space one can also 
introduce operators, which will be called superoperators. 

The time evolution of the spin density matrix p^ sat- 
isfies the Schrodinger equation (hereafter we put = 1, 
i.e. use frequency units for energy parameters) 

Pi, = -iH^pi,. (3.3) 

In this equation the spin Hamiltonian superoperator 
(operator in the Liouville space) expressed in terms of the 
Hamiltonian Hi, (in the Hilbert space) by the relation 

H„py = Hu\pu) Hi^pi, - p„H„. (3.4) 

For matrix elements of Hi, we get {i'p,^^'^\H,^\h'^^^'^) = 
(j/^Ji7^|j/^JViMi - {i^f,'jH,,Wf,'JS^,^f,^. 

IV. SPIN HAMILTONIAN AND SPIN 
RELAXATION 

In general the spin Hamiltonian iJ^(r^), (i' = e,h), 
depends on the coordinate r^, of i^-polaron, due to in- 
homogeneity of the medium. This dependence results 
in fluctuations of H,^{t) = H,^(ru{t)) caused by hopping 
motion of j/-polaron, i.e. stochastic changing rj,(i). The 
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Hamiltonian can conventionally be represented as a sum 
of the stationary and fluctuating parts: 

H,{r,{t)) = H,^ + H,^ {v,Xt)) + H,^ (r,(t)). (4.1) 

in which the stationary part 

H,^^ = .g„/3(S, • B) = with = g^^B (4.2) 

is the average Zeeman interaction of the electron spin 
with the magnetic field B. The fluctuating parts 
H^^{r^{t)) and H^^{r^{t)) represent HFI and AZI con- 
tributions, respectively. 

Fluctuating HFI and AZI induce spin relaxation, 
which can be treated within the short correlation time 
approximation!^ In this approximation the relaxation is 
described by the Bloch-Redfield equations for spin den- 
sity matrices of polarons.~ The relaxation kinetics is de- 
termined by the relaxation supermatrices of the form 

Wq^=wP^PP + wIP:\ {Q^H,A), (4.3) 

for both HFI (H) and AZI {A) mechanisms, where 

P}^ = (|z.,J-|z.__))((^^J-(i._J), (4.4) 
C = W,_){iy^J + W_^){iy_,\. (4.5) 

are the operators in the subspaces of diagonal (PP) and 
non-diagonal (P") elements of the density matrix. The 
first and second terms in eq. (|4.3p describe the popula- 
tion and phase relaxation, respectively, with rates^ 

%'„=%.'^«(^'')' (4-6) 
Here is the numerical parameter, depending on the 
relaxation mechanism (specified by Q — H, A), iD^ is 
the characteristic rate (see below), and 

/•oo 

(^) = ^qT' / % (*) (4-7) 

with 

/■oo 

- / dt%it) (4.8) 
Jo 

is the normalized Fourier transformed correlation func- 
tion satisfying the relation Ji/(0) = 1. 

The form of functions J^i^j) depends on the mecha- 
nism of fluctuations. In the considered model of polaron- 
hopping induced fluctuations these functions are essen- 
tially determined by hopping kinetics and specific fea- 
tures of the orientational distribution of molecules re- 
sponsible for the HFI and AZI. In our work we will need 
only the most general properties of J^((jj). They can 
be understood in a simple model of random Hamiltoni- 
ans H^^{r^) and H^^{r^) uncorrelated at different sites 
(with zero mean values), in which ^^it) ^ = 
(exp(-w;^^^Jt|))^^^^ and 

^:(.) = (^.^.^ «,^. +u.rX^^ /«.^. l^^^ . (4.9) 



A. Relaxation mechanisms 

1. HFI and HFI induced relaxation 

The HFI is determined by the spin-spin interac- 
tion of the electron and paramagnetic nuclei Vj (with spin 
/i,^. ), localized in close surrounding of electronic spins Se 
and Sfi of e- and /i-polarons. In the realistic case of a 
large number of nuclei the interaction can quite accu- 
rately be approximated by that of spins Se and Sh with 
(classical) random magnetic fields Be and B/j, respec- 
tively, whose distributions are isotropic and Gaussian 
with mean squares {B"^) ^ I^. (1 -I- 1^. )al. , = e,h), 
determined by hyperfine coupling constants a^,^ 

Hopping of e- and ft,-polarons results in sudden chang- 
ing of nuclear magnetic field, which can be considered 
as a stochastic vector B^(i), = e,h), with (B^) = 
and the correlation function of projections Bi, (t), (g = 
X, y, z): {B,^ {t)B,^, (0)) = \5,,. {Bl)^l (t). 

The HFI mechanism predicts the relaxation superop- 
erator Wj^ of the form (|4.3I) with the rates and 
given by eq. (j4.6l) . in which 

= 1 and w^^ = i(5./3)'(S^)T„„, (4.10) 

with the correlation time defined in eq. (|4.8p . 



2. AZI and AZI induced relaxation 

The AZI results from the deviation of ^jy-factors 
of e- and /i-polaron spins from the free electron value 

= 2j^i^ In general, g-factors are actually ^jx-tensors — 
whose eigenvectors are determined by the geometry 
of molecules, at which e- and ft.-polarons are located. 
The gi/-tensors are conveniently represented as sums of 
isotropic {gv) and anisotropic {5gu) parts: 

9y = 9u + g'^, where = Tr(g,,). (4.11) 

The representation (|4.1ip implies that Tr(g^) = 0. Usu- 
ally, for systems under study eigenvalues of gj^-tensors 
are close to go = 2: 6gy = g^ - g^ ^ \\g'^\\ < 10"^. 

The AZI part of the interaction of the electron spin v 
{v = e, h) can be defined as 

Ha_^ = /3S,5,B - Hz^ = pS.glB. (4.12) 

In the AZI induced relaxation mechanism the rate su- 
peroperator VF^ is written as (j4.3p with rates and 
w" (|4?6l) . in whicbSi 

p^ = I and w^^ = ^{gl: gl){l3B)\^^, (4.13) 

with T^^^ given by eq. (j4.8|) . 
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B. Hamiltonian and relaxation in polaron pairs 

The PP spin Hamiltonian is, in general, a complicated 
function of coordinates Vh and Fg, determined by the elec- 
tron spin dependent PP interactions: spin exchange and 
dipole-dipole interactions. Fortunately, the their effect 
on majority of MFEs is weak and can be neglected. In 
such a case the PP spin Hamiltonian can be represented 
as a sum of spin Hamiltonians of e- and /i-polarons: 

H = H,^+ H,^ = -gPBS, + ^Agf3B{S,- S„J, (4.14) 

where = S'g + S"^ , 

9^\{9e+9h) and Ag^g^-g,^. (4.15) 

Within similar assumptions, the PP spin relaxation 
can be described by the superoperator, which is also a 
sum of corresponding operators for separate polarons: 

W = We + Wh^y] (W^ (4-16) 

The PP spin evolution is described in the four-state 
Hilbert space. In principle, one can use any basis in this 
space, for example the basis of states |ep)|/ip'), {fi, fi' = 
±), of the pair of non-interacting polarons. However, in 
what follows it will be more convenient to use the basis 
of eigenstates of the z-projection, S^, of the vector S: 
IS"), |r^), {fi = 0,±), which correspond to the total; spin 
5 = (5) and 5 = 1 (T): 

1^) = j^i\e+)\h-)-K)\h+)), (4.17) 
l^o) - ^(|e+)|/i_) + |e_)|/.+)), (4.18) 
= \e±)\h±). (4.19) 

In addition to the states in the Hilbert space it is worth 
introducing the states in the Liouville space 

\XY)^\X){Y\, iX,Y^S,T^). (4.20) 

V. STOCHASTIC LIOUVILLE EQUATION 

The specific features of PP recombination processes 
(both geminate and bulk) can, in general, be expressed 
in terms the PP density matrix p{r, t) . In the case of 
Markovian e—h relative motion p{v, t) is known to satisfy 
the SLE,^^'^^ which within the diffusion approximation 
is written as 

-[k + k{r)+L]p, where k = iH + W, (5.1) 

H = , . . . ] is the superoperator representation of the 
PP spin Hamiltonian H [see eq. (j4.14p ] in the Liouville 
space, and W is the superoperator of PP spin relaxation, 
defined in eq. (14.161) . 
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FIG. 1: Schematic picture of the PP interaction potential 
u(t) = U{r)/ (kBT), in which d is the distance of closest ap- 
proach, ri, is the coordinate of the well bottom, Ui, = —u{rb) 
is the well depth, and 1^ is the Onsager radius. 

The relative diffusive motion of e- and ft,-polarons is 
described by the Smoluchowski operator L, defined by 
the expression 

L/9 = -L'Vr(VrP + pVrU), (5.2) 

in which r = — the relative PPq coordinate, 
D = Dh + De is the relative diffusion coefficient [see 
eq. (|2.ip ]. and u{r) = U{r)/ (kBT) is the dimensionless 
PP interaction potential assumed to depend only on the 
distance r = |r|. 

In our analysis we will consider the spherically sym- 
metric problem, thus reducing it to studying the evolu- 
tion along the radial coordinate r = |r| = |r^ — r^l only. 

The potential u{r) is assumed to be of the shape of 
potential well (see Fig. 1). This potential is character- 
ized by the barrier, of height u,. = u{d) — u{r\,) > 1, 
at the distance d of closest approach of polarons (which 
models the suggested smallness of the reaction rate) , the 
coordinate ri, of the bottom {n, > d), the well depth 
Uj, = —u(ri,) > 1, as well as the Onsager radius 

/•OO 

lc= drr~^e''^'-\ (5.3) 

The term K{r)p describes spin selective recombination 
assumed to occur only in S'-state of the pp>i2i2£ 

k(r)\p) = ln{r){Psp + pPs) ^ n(r)Vs\p), (5-4) 

where K{r) is the distance dependent recombination rate. 
In eq. (|5.4I) Ps — |S')(S'| is the operator of projection on 
the S'-state and 

is the superoperator, controlling the spin dependence of 
reactivity (and accompanying dephasing), in which 

Pxy = \XY){XY\, {X,Y^S,T^), (5.6) 

are superoperators of projection onto the states (|4.20p 
in the Liouville space. The dependence K(r) is typically 
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short range and does not affect MFEs very much, so that, 
for simphcity, one can apply the contact reactivity model 

K{r) ^ Kj{r-d), i.e. k{r) ^ KJ{r~d) (5.7) 

with K^=n^V^. 

The polarons are assumed to reflect at the distance 
of closest approach d. The process is described by the 
reflective inner boundary condition (VrP+pVru)^^^^ = 0. 
In the model of contact reactivity (|5.7|) . however, one can 
formally reduce the effect of reactivity to the modification 
of the inner boundary condition: 

[WrP + pWrU~{KjD)p)]^^^=0. (5.8) 

As for the outer boundary condition it is different for 
geminate and bulk processes (see below). 

VI. GEMINATE PP RECOMBINATION 

The kinetics of geminate PP recombination will be 
analyzed assuming that PPs are created with isotropic 
distribution, localized at a distance r — Vi, and in the 
spin state, determined by the density matrix Pi{r). This 
means that the SLE (|5.ip should be solved with the ini- 
tial condition 

p,{r) = p(r, i = 0) = (4^r2)- V/(r - r,). (6.1) 

The outer boundary condition at r — ^ cx), is written as 
p{r oo,t) 0. 

Spherically symmetry of PP interactions, as well as 
spherical symmetry of initial and boundary conditions 
ensures that the spin density matrix is also spherically 
symmetric: p{r,t) = p{r,t). In this case the PP recom- 
bination yield Yr, can be expressed by formula 

Yr = j dt Jd^rTr[K{r)p{r,t)] 

= 4TrKsd^S\p„{d)\S) = ATTKsd^SSlpoid)), (6.2) 

where the matrix pQ{r) is the Laplace transform Pe{r) = 
/(J dte~'^*p{r,t), evaluated at e = 0: Poir) = p^^air). 

The matrix p^ir) satisfies the steady state variant of 
the SLE (EU 

{A + Lr)po^ p,{r), where A = iH + W , (6.3) 

and Lr is the radial part of the operator L defined by 

Lrp = Dr-^Vr[r^{VrP + pV^it)]. (6.4) 

Noteworthy is that in the SLE (|6.3p we have omitted 
the reactivity operator K{r) (|5.7p . expressing its effect 
by the properly chosen inner boundary condition (j5.8p . 

For the initial condition (|6.ip the solution of this equa- 
tion is easily expressed in terms of the Green's function 
of the equation (|6.3I) 

Go = (A-l-L^)-i : (6.5) 



Yr = >,s{d/nf{SS\Go{d,n)\pg). (6.6) 

In general, G'o(r, r^) can hardly be obtained analyti- 
cally. However, fairly simple analytical expressions for 
GQ{r,ri) (and thus for Yr) can be found in some impor- 
tant particular cases. 

A. Freely diffusing polarons. 

1 . General formulas 

In the absence of PP interaction potential [u{r) — 0] 
the solution of eq. (16.51) can be obtained in a matrix 
formi^i^ For simplicity we will assume that PPs are 
created at a distance of closest approach, i.e. — d. 

The analytical expression can conveniently be obtained 
with the use of the representation 

G„^D-\r,/r)g^ (6.7) 

In this relation 

9o = (fco - V^)-\ with fco = {A/Dy^\ (6.8) 

is the one dimensional Green's function, which satisfies 
the inner boundary condition similar to ([5.8^ : 

(V, - qs)9Q\r=d = with Qs = + (KJD). (6.9) 

The function is given by formula^^*^ 

5o(r,r') = [e-''^oir-d)X + e-f^o\r-r'\^(^2k,)-\ (6.10) 
where 

A = (^ + ly^e ~ l)e-''o('-'-'^^ and § = q-'^ka. (6.11) 

Substitution of this formula into eqs. (j6.8|) and (|6.2p 
results in the expression for the recombination yield: 

Y^Yf^ TviPsiPfPg)] ^ {SS\ P^lpg), (6.12) 

in which 

Pf=rr,il+dkorej)-\ (6.13) 

is the spin dependent supermatrix of reaction/relaxation 
probabilities. In eq. (I6.13|) 

Vr^ = i/d and Ve^ = I - l/d, (6.14) 

are the supermatrices of probabilities of reaction and 
escaping (respectively) of PPs, created at = d, in 
the absence of PP spin evolution. The probabilities 
are essentially determined by the supermatrix / of re- 
action/relaxation radii, corresponding to the rate super- 
matrix i4r(r) [given by eq. (j5.7[) ]: 

i^d-q-^ = djf/{l+jf) with = dkjD. (6.15) 

Formulas (|6.12p - (|6.15p for the yield of recombination of 
freely diffusing polarons are seen to reduce the problem 
of evaluation of MFEs to simple matrix operations. 
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2. Fast P-diffusion limit 

The limit of fast relative diffusion (or slow spin evolu- 
tion), when c?||fco|| ^ 1 , — is of special interest for our 
further analysis. In this limit can be found in the 
approximation linear in dka: Pj: « ~ dk^Pef), i-C. 

Yf « {lssld)p, ~ lss{SS\hV,^\Pg). (6.16) 

where Iss is the reaction radius in S'-state and — 
{SS\pg) is the initial population of this state. 

For the model (|4.14p . (|4.16|) the yield Yf can be ob- 
tained for any initial state \pg) and, in particular, for 

\Pg) ^Ps\Ps) +Pt\Pt)^ (6-17) 

in which and p^ are the probabilities of population of 
S and T states, respectively (pg + P^. = 1); 

\p,) = \SS) and \p,) = i EJ^a^^m) (6-18) 

Calculation of Yf reduces to evaluating the superma- 
trixA:o = {iV^Y^ dtt-^/"^ {l-e'^^^) . Substituting 
thus obtained /cq into eq. (|6.16p . one gets 

Y « ^ [Ps + T2(*J(P.^; - 3Ps^4)] , (6-19) 
where "P^^ = 1 - lj,j,/d =\,V^^. = l- P^. = Z^^/d, and 

= + 2Re(^w„-HiAw )] /\/Z). (6.20) 

In this formula 

Ac^ = (ge - gh)liB = Jg/3B. (6.21) 

is the difference of Zeeman polaron frequencies and 

Wj, = 2[wl + < ) and u.„ = u;« + < , (6.22) 

are the total rates of population (wp) and phase (w„) re- 
laxation in the PP, which are the sums of contributions of 
the HFI and AZI induced relaxation rates. These contri- 
butions are, in turn, the sums of corresponding relaxation 
rates in e- and /i-polarons: 

= "^H. + "^H,, ; ^ "^le + "'Ift ' (9 = -P; (6.23) 

The total rates and w„ are represented as sums of 
HFI and AZI induced relaxation rates to clearly reveal 
two contributions with different i?-dependence: decreas- 
ing [w^ {B)] and increasing [w^ {B)] as B increases. 

B. PP recombination in the presence of interaction 

The PP interaction potential u{r) can essentially af- 
fect the MFEs on PP recombination. For pure repulsive 
interaction [u{r) > 0] no significant effect is expected 



except some change of reaction and relaxation radii, i.e. 
the elements of the supermatrix I introduced in eq. (|6.15p 
for freely diffusing polarons. Much stronger effect is pre- 
dicted in the case of attractive interaction, especially for 
potentials of the shape of potential well (Fig. l)j^ In 
what follows we will discuss the MFEs in the special case 
of deep potential well with 3> 1. 

The analysis can be made with the use of the recently 
proposed method of rigorous analytical analysis of the 
kinetics of diffusion in the presence of the wellj^S, One of 
important results of this analysis consists in the fact that 
in the limit of deep well (corresponding criterion is given 
below) the exact SLE (15. ip is equivalent to the model of 
two kinetically coupled states: the state within the well 
(the diffusive cage state), at d < r < Ic, and the state 
of free diffusion outside the well, i.e. at r > l^^^l The 
spin/space evolution in these two states are described by 
density matrices 

n{t) ^ An drr'^p{r,t) and a{r,t) = rp{r,t), (6.24) 

Jd 

respectively, satisfying equations^ 

n = [SYK+<y{lc,t) - {K- + Wr)nl (6.25) 
cr = [DVla + {SiK^n - K+a)6{r - Q], (6.26) 

in which Si = (47r/c)~^. The terms proportional to 
K± represent the kinetic coupling (transitions) between 
the two states, with transition rates K± satisfying the 
relations)^ K± — )• 00 and K+/K^ = Kg = Z^^, 
where 

= [ drr^e-''^'^^ (6.27) 

is the partition function for the well. 

Equations (|6.25p and (|6.25p will be solved assuming 
that the PP is initially created within the well, i.e. 

n(0) = pg and cr(r, 0) = 0. (6.28) 

As to the boundary conditions for c{r,t), they are given 
by Ic'^rO'ir, t) — cr(r, i)|r=/e — and a{r — >■ 00) = 0. 

The term WrU in eq. (|6.25p describes the first order 
reaction in the well with the supermatrix of rates 

and % = {dk J D) {dK e -'''^'^^ ) 

Solution of eqs. (|6.25p and (|6.26p by the Laplace trans- 
formation in time leads to formula for the recombination 
yieldi'* similar to cq. (I6.12p : 

Y, - T,[Ps{P,Pg)] ^ {SS\ PJp,), (6.30) 

but with the supermatrix P^- replaced by 

P, = Wr[W, + k + We{lcko)Y^, (6.31) 
= VrA^ + {k/W,+lchYeX\ (6.32) 
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in which K = iH + W, fco = (A/D)i/2, and 

Wc^Wr+ We, with We = Dle/Z^, (6.33) 

is the supermatrix of total cage decay rates, represented 
as a sum of the rate Wr of reaction in the weU (cage) and 
the rate We of escaping from the weU.— These rates 
determine the probabihties of reaction in the weh iVr^) 
and escape from the weU [Ve^ ) : 

Vr, = Wr/We and Ve, = We/We- (6.34) 

The detailed analysis showa^l' that expressions (j6.3ip 
and (|6.32p are valid in the limit of not very fast PP spin 
evolution, or not very large size of the well Sc = Ic — d, 
when (5c||A:o|| ^ 1 (this inequality does not, in principle, 
mean that ^dl^oll !)• 

According to obtained formulas, the effect of the well 
manifests itself in the formation of the diffusive cage, 
whose evolution is described as the first order reaction 
(with the rate Wr) and escaping (with the rate We). 
This simple first order kinetics is, however, perturbed 
by the contribution of particles escaped but recaptured 
back into the well, resulting in the term ~ (Icko) in eqs. 
(|6.3ip and (|6.32p . which gives negligibly small contribu- 
tion to MFEs for deep wells (when We is so small that 

|1A|1/We > ^cllfcoll), but strongly affects MFEs in the 
opposite limit of large escaping rate We- Moreover in 
the limit of fast cage decay (or slow PP spin evolution), 
||A||/PVe ^ /c||fco||, the obtained formula reduces to that 
dO'^ for free diffusion,^ 

VII. BULK PP RECOMBINATION 

The MFEs in bulk recombination are also described by 
the SLE (15.11) but with the boundary condition at r — > oo 

pir ^ oo,t) = ^ (7-1) 

corresponding to the homogeneous spatial distribution of 
and the equilibrium spin states of polarons [represented 
by the unity matrix = E;.=s,To,± 1/^) (/^H- 

In the case of bulk recombination the observable under 
study is i?-dependent recombination rate K(t). For sim- 
plicity we will discuss the static value K = K{t oo) 
which can be represented in the form very similar to that 
of the expression (|6.2|) for the recombination yield: 

K = J(fr TT[k{r)p{r, t oo)] 

= ^nKs(f{S\po(d)\S)=ATTK,d''{SS\Md)), (7.2) 

This formula is written, taking into account that the 
equation for the static solution Pstir) coincides with the 
Laplace transform Pe=o(^) = Po('") [however, with the 
outer boundary condition (|7.ip ]: Pstif) = Poi''')- 



A. Freely diffusing polarons. 

The stationary solution Po(^) of the SLE (|5.ip satisfy- 
ing eq. ()6.3p with pi — can be found by solving more 
simple equation for crQ(r) — r~^pQ{r) 

(kl - V^)5o = with fco = (A/D)^^^. (7.3) 

The function CTq (r) satisfies the inner boundary condition 
(Vr — qs)(^o\r=d = Oj similar to (|6.9p . and the outer one, 
which, according to eq. ()7.ip . is written as 

ag{r ^ oo) ^ rp^. (7.4) 

By direct substitution into eq. ()7.3p one can show that 
the solution of this equation is given by 

ao{r)^[-e-^'^^'^Hl + 0)-H+r]p,. (7.5) 

where 6 is defined in eq. (16. lip . With the use of this 
solution the recombination rate is represented as 

K^Kf^ K'j-Yf^, where = {SS\ Pf\p,}, (7.6) 

In this formula Yf is the geminate recombination yield 

for the equilibrium initial Pg — p^, P/ is the superop- 
erator of reaction/relaxation probabilities defined in eq. 
(|6J3)) . and K'j = AnDd is the rate of PP contacts. 

Note that the dimensionless rate Kj:/K^ can be related 
to the yield Yf^ of the geminate PP recombination for 
(triplet) initial density matrix pg = p^.'. 

Kf/K'I^Yf^^lVr'^il + SYf^), (7.7) 

where V^^ = Igs/d is the PP recombination probability 
in S state. Equation ()7.7p is obtained using the relation 

= Vrj{l + dk^Ve^)-^ = Vr^ - P/dA^Q^e^. 

B. Interacting polarons 

In the considered case of attractive interaction the ex- 
pression for the PP recombination rate can be ob- 
tained within the two state model (Sec. VLB). For bulk 
reactions the spin density matrix, represented in terms 
of spin density matrices of the state within the well Uq 
and the free diffusion state CTq (r) , satisfy the steady state 
variant of eqs. (|6.25p and (|6.26p . i.e. equations with 
'^0 = '''o = 0; but with the outer boundary condition: 
cro(r -> oo) = rp^. 

Solution of these steady state equations leads to the 
following expression for the PP recombination rate 

K^K, = K^e, with Ye^ = {SS\ PJpJ, (7.8) 

in which Ye^ is PP recombination yield for the equilib- 
rium (spin) initial state, P^ is given by eq. (|6.3ip . and 
= AttDIc is the rate of capture into the well (cage). 
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The relation between the dimensionless rate K^/K^ 
and the recombination yield [for the triplet (T) initial 
condition], similar to eq. ()7.7p . can also be derived in the 
case of attractive interaction: 

KJK', = Y,^ = \V^^ (1 + J , (7.9) 

where V^:^ = {SS\{WrlWc)\SS) is the probability of in- 
cage reaction in S state. The derivation is based on equa- 
tion = Vr, - Pj\ k/We + {lcko)]Ve,, with Vr, and Ve, 
defined in eq. ()6.34p . 

VIII. RESULTS AND DISCUSSION 
A. General remarks 

1. Validity of approaches 

Before the analysis of the MFEs it is worth adding 
some comments on validity and accuracy of the proposed 
mechanisms of spin relaxation. In the applied model the 
relaxation, assumed to be induced by fluctuating HFI and 
AZI, is described by simple Bloch-type equations, which 
are valid at times longer than the correlation times 
and Ta of HFI and AZI fluctuations. This means that the 
relaxation mechanisms are applicable only if the time of 
MFE formation ^ Th, Ta- 
in the free diffusion model the MFE-formation time 
Tg is completely determined by the time of spin evolu- 
tion in polarons: r_g ~ ||^||~"^ ^ ''^n^ ^^p^ A'^^)^^ 
Taking into account that the relaxation rates (w^ and 
w^) and coherent evolution frequencies {'^ Aw)] are typ- 
ically of order of (or less than) lO^s"^?^^ one obtains: 
Ts > 10"* s. 

Similar estimation for is valid in the diffusive cage 
model. The additional time parameter, which could 
change it, is the characteristic inverse rate of the cage 
decay Tc = llW^dl^^- However, this time is also fairly 
long: Tc > i/7(^^e"", where wq is the characteristic hop- 
ping rate and Ua is the activation energy for escaping 
and/or reaction processes. The assumption of deep well 
{ua ^ 1) leads to the estimation Tc ^ w^^. 

As for correlation times they are expected to be of or- 
der of inverse hopping rates: t^ ^ Ta ^ '"'0'^' since HFI 
and AZI fluctuations result from stochastic polaron hop- 
ping. Typically the rate Wq ~ 10^ — 10"'^^ s~^, which cor- 
responds to room-temperature mobilities fie,h ~ 10~® — 
10'^cmV(Vs),JI therefore we get Th,a 10"^ - 10-"s. 

Comparison of these estimations shows that the above- 
mentioned validity criterion, t^^Tc ^ Th,Ta, is fulfilled 
for a large number of semiconducting systems considered. 

Concluding the discussion note that the validity of the 
diffusion approximation implies negligibly small contri- 
bution (to the MFEs) of small times t ~ tq = Wq^, at 
which the diffusive motion is not yet formed. In the con- 
sidered /- and c-models, according to the above relations, 
their MFE contributions SY,, = Y^{B) - r^(0), {p = 



f,c), can be written as SYf ~ {wQTg)~^^^ and SYc ~ 
t^/Tjj ^ e"''/(wgrg). As to the contribution SY^^^ of 
small times t ^ tq, it is, evidently, represented by 
JYgj ~ (wgTg)"^, i.e. in the considered limit WqT^ <^ 1 
the contribution SY^^ is really small: SY^f <C SYj:, SY^ and 
the diffusion approximation is applicable. 

2. Parameters and observables 

In the proposed models B-dependence of the recombi- 
nation yield Y(B) is determined by a very large number 
of parameters of the model. To reduce this number in 
our illustrative discussion of most important properties 
of Y{B), we consider the particular (but representative) 
variant of the model, in which population and phase re- 
laxation rates {w^ and w", respectively) are the same in 
e- and ft,-polarons: 

with Q = H,A; and 

<=%^qH; <=PqK+K- (8-2) 

In formulas (|8.ip and (18. 2p the parameter Q denotes the 
contributions of HFI (Q = H) and AZI (Q = A) mech- 
anisms, and the parameters and are defined in 
eqs. ()4.10|) and (|4.13p . In these formulas we neglect the 
difference of Zeeman frequencies in functions ^/^(a;^), i.e. 

took UJe ~ UJfi ^ ~ fjPB. 

The MFEs on geminate and bulk PP recombination 
are found to be closely related [see eqs. (j7.6|) - (|7.9|) ]. so 
that it is sufhcient to analyze, for example, the MFE 
on bulk process, i.e. i3-dependent (dimensionless) rate 
K^{B)/K^ — Yf_i^{B). In what follows we will discuss 
the function 

y^iB) = [Y^jB)^Y^jO)]/Y^jO), (m-/,c). (8.3) 

In our work we have obtained fairly simple matrix ex- 
pressions for the MFE y^ (B) in two models of relative 
polaron motion. For qualitative understanding of the 
properties of these functions, however, it is of certain in- 
terest to get simple approximate analytical expressions. 
Below we will derived them in some limiting models. 

B. Simple limiting models 

1. Fast freely diffusing polarons. 

In the case of freely diffusing polarons of special inter- 
est is the limit of fast polaron diffusion or slow polaron 
spin evolution, in which |(fcolM ^ 1- 

General consideration of the problem predicts for the 
recombination yield the expression (|6.19p . In what fol- 
lows we will mainly discuss the shape of the MFE Yf^ , 
which can be written as 

« 3^^, + ll(nV'(rffc»), (8.4) 
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where fc^ = [y/W^ + 2Re(-y/w„ + iAuj)]/VD [see eq. 
(|6.20p ]. V^^ = Iss/d is recombination probability in S 

state, Wp and Wn are population and phase relaxation 
rates defined in eqs. (|6.22p and (|6.23p . respectively, and 
Aw = {ge — gh)PB = {Ag)[3B is the difference of Zeeman 
frequencies [eq. ()6.2ip ]. According to formulas (j6.22|) 
and (j6.23|) the rates Wp and Wn are the sums of the cor- 
responding rates for e and h polarons, which are, in turn, 
the sums of contributions of HFI and AZI mechanisms: 
Wp = 2{w^ + w^) and u;„ = + so that 



WpiB) 



= 2[w^J^{B)+wJ,B)JJ,B)], (8.5) 
+ + \wp{B). (8.6) 



w. 



In these equations we have introduced the function wJ^B) 
to emphasize that (unlike ) lu^ depends on B [see eq. 
Km ]: wjB)^B\ 

Noteworthy is that, although formula (18.41) is derived 
in the fast diffusion limit dkw 1, it appears to be quite 
accurate (within 15 — 20%) even at dk^ ^ 1. 



2. The limit of weak reactivity in the well. 

Majority of specific features of MFEs in the presence of 
the well (diffusive cage) can be analyzed with the limit 
of weak reactivity, when Q = ||M^r||/M^e ^1- In this 
limit the expression for the yield Yc can be derived by 
expansion of the general formulas (|F (lOTD . and (USD 
in small Cc- In the lowest order we get 



i-ps 



+^{Vlf{N{wp) + 2Re[N{w,,+iALo)]}, (8.7) 



where = {SS\Wr/Wc\SS) < 1 is the probability of 
reaction in the well and 



N{e) 



1 



with ^ = l^yWjD. 

It is seen that in the limit of slow spin evolution, when 
Wp, |w„ + iAuj\ <C £,^We, we get N{e) ^ ^/e and formula 
(18. 7p reduces to eq. (|8.4p . derived in the free diffusion 
model. Of great interest is also the opposite limit, in 
which the evolution kinetics is described by simple first 
order kinetic equations and iV(e) ~ (e/Wg)[l-|-(e/lVg)]~^. 
This kinetics shows itself in the typical analytical depen- 
dence of MFEs on relaxation rates and the splitting Aw. 
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FIG. 2: Dependences y^(&//) on = B/B^ for the HFI 
induced relaxation mechanism, calculated in free diffusion 
{fi — f) and cage (/i = c) models: for Lorenzian (Figs 
2a and 2b) and non-Lorenzian (Fig. 2c) shapes of J7^(_B). 
Fig. 2a shows yf{bu) for different 77^ = d^w^ / D: 77^ — 
0.1 (1); 0.3 (2); 1.0 (3); 3.0 (4). Fig. 2b shows yAbn) for 
Tc = w^/Wc = 0.5 and different = Wc/Wu'. = 
5.0 (1); 1.0 (2); 0.1 (3); 0.01 (4). Fig. 2c shows arbitrar- 
ily scaled y^ibg) for non-Lorenzian J^[B) (Sec. VIII.C.l) 
with ri^ = 10"'' for = / (1) and r„ = 1.0, = 5.0 
for /i = c (2). Circles represent the prediction of formula 
y^^{B) ~ [S/(Bo + |S|)]^ with Bo = 2.35B„, obtained by 
fitting y^ibn). 



Here is defined in eq. (|5.5I) and w^. = D/{Kj.Z^). 

We also discuss in detail the special case of Lorenzian 
shape of functions Jq{B) = (1 -f- S^/S^)-!, for which 
the relaxation rates are conveniently represented as 

< = wj{l + hi) and < = w*yj{l + 62 ), (8.10) 

where — B/B^ with B^ being the width of J [B). 



In our analysis (to reduce the number of parameters) 
we consider the limit of high reactivity at a contact, 
K^/D 3> 1, in which reaction/relaxation supermatrices 
I and are independent of : 

i « dV^ , W, « w,V^ with = 2V^ - 4, . (8.9) 



1. HFI induced relaxation 

Specific features of dependences yp[B) for the HFI- 
induced relaxation mechanism are demonstrated in Fig. 
2 for the free diffusion (/i = /) and diffusive cage (/i = c) 
models. In both models this mechanism leads to the de- 
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creasing functions y^{B). The decrease results, evidently, 
from that of population relaxation rates as B increased. 

As expected from analytical expressions (|8.4p and (18. 7p 
these models predict similar dependence y^{B). In a wide 
region of values of the diffusion coefficient and reactivity 
the change of these parameters is found to lead mainly 
to the change of the MFE amplitude rather than to that 
of the shape of y^(i?)-dependence. 

The shape is mostly determined by the mechanism of 
relaxation, or more accurately by the form of the cor- 
relation function J^{B). For example, in the case of 
Lorenzian (B) [see eq. (|8.10p ] the MFE shape is close 
to Lorenzian with the width of order of that for J'^ {B). 

In our short discussion we are not going to review all 
properties of the shape of yAB). Majority of them can be 
described with formulas (|8.4p and (|8.7p . More thoroughly 
we will only analyze the possibility of interpretation of 
the non-Lorenzian MFE shape y^xiB) ^ [B /{B^ + \B\)f , 
observed in some experimentsi^ii^ 

Note that for the non-Lorenzian function J^{B) the 
MFEs y^{B) are non-Lorenzian as well. Just this prop- 
erty of the HFI induced relaxation mechanism allows one 
to describe the observed dependence y^^{B). The re- 
action yields yf{B) and y^{B), shown in in Fig. 2c, are 
obtained for J^{B), — f,c), represented as sums 
of Lorenzian contributions Lo'^{bjj) = [1 + ^^/(jn)^]~^, 
in which j„ are the numerical coefficients: J^^ (B) 

J2n=i '^ra^On(^/f )) = /; c), whcrc af are the weights of 
the Lorenzian contributions, taken the same both for the 
free diffusion and diffusive cage models: a{ = a'l — 1, {i — 
1, . . . ,6). The values of j( and jf, conveniently repre- 
sented as vectors j'^ = (jf , ■ . ■ are different in these 
models: = (1, 2, 4, 4, 6, 10) and j'^ = (1, 3, 4, 5, 7.5, 20). 
The chosen number of terms, n = 6, seems to be the 
smallest of those, which are sufficient to get the yields 
quite close to yg^{B). 

Thus calculated yf{B) and y^B) agree with y^xiB) 
quite well in the wide region of B values. 

The obtained wide spectrum of widths, required in 
both models, is qualitatively consistent with the wide dis- 
tribution of hopping rates (quite natural for disordered 
semiconductors), which determines the behavior of cor- 
relation functions, according to eq. ()4.9|) . 



2. AZI induced relaxation 

Unlike the HFI induced relaxation the AZI induced 
one results in the increasing dependences y^{B) in both 
models of PP migration. As for the characteristic fea- 
tures of the shape of these dependences, they are not 
very similar to those for HFI induced relaxation. Typi- 
cal dependences yf{B) and yc{B) are displayed in Figs. 
3a and 3b for Lorenzian shape of Jj^{B) [see eq. (|8.10|) ]. 

It is seen from Figs. 3a and 3b that for AZI mechanism 
the widths of yf{B) and ydB) are larger than those for 
HFI mechanism and, in general, are independent of the 
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FIG. 3: The dependences y^ib^) on = B/B^ for the AZI 
induced relaxation mechanism and for the Lorenzian shape of 
J^{B), calculated in free diffusion (/i = /) and cage (/i = c) 
models (Figs. 3a and 3b), and j/c(^ah) for the HFI+AZI 
relaxation mechanism [assuming J^[B) = J^{B) and 6^ = 
bfj;] (Fig. 3c). Fig. 3a displays y;{bj<^) for different -q^ = 
d^f^lJD: = 0.02 (1); 0.2 (2); 2.0 (3); 10.0 (4). Fig. 3b 
displays ydbn:) for = w^/W^ = 0.5 and different = 
Wc/w*^: K,^ = 0.05 (1); 0.5 (2); 5 (3); 50 (4). Fig. 3c displays 
yc{bjy) for = 1, ~ ^c/w^ — 1, and different = 
Wc/w*^: = 1 (1); 2.25 (2); 4 (3); 10^ (4); 10^ (5). 



width B^ of J^{B). The reason of this independence 
consists in that, unlike the HFI induced relaxation, the 
AZI induced one results in dephasing, whose rate rapidly 
increases with B: ^ w^^ ^ B^ [see eq. (|8.2p ]. In 
such a case the width of y^{B) is mainly determined by 
the kinetic saturation of the MFEs observed for large 
dephasing rates [Sec. IV. A2 and eq. (|8.ip ]. In the 

free diffusion model the saturation occurs at \k^\d > 1, 
i.e. outside the region of validity of the fast diffusion limit 
(|8.4p (Sec. VI. A2). As to the cage model, the saturation 
in this model is observed for [|A[|/[|W^[| > 1. Within 
the weak reactivity limit (18. 7p this condition reduces to 
simple inequality w„/W^ > 1. 

Most clearly the peculiarities of MFEs, resulting from 
AZI induced relaxation, manifest themselves in the lim- 
its of fast free diffusion and fast escaping from the cage 
(wp^/W^ <C 1). In these limits obtained formulas pre- 
dict [see eqs. (|8.4p and (18. 7p ] monotonically increasing 
behavior: yf{B) ^ \B\ and y^iB) ~ B^ in the region of 
B wider than that {B < B^) expected from eq. (|8.10p . 

It is worth noting that the dependences close to 
y{B) ~ \B\ has been recently observed in a number of 
experiments ji^ii^i^ The proposed analysis allow for quite 
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reasonable interpretation of this, at first sight, strange 
behavior. 



3. Superposition of HFl and AZI induced relaxation 

In the case of comparable contribution of HFI and 
AZI induced relaxation the free diffusion and diffusive 
cage models predict a large variety of MFE dependences 
y^{B), (/i = /, c). In our brief analysis it is hardly possi- 
ble to describe all types behavior. In general, the 
superposition of HFI and AZI (HFI-I-AZI) induced re- 
laxation mechanisms results in non-monotonic functions 
y^{B). Some representative examples of them, as applied 
to VciB), are displayed in Fig. 3c. Analysis shows that 
the resulting shape of y^{B) dependence is essentially de- 
termined by the relative widths and relative weights (i.e. 
= w^/Wc and = zD* /Wc) of both mechanisms. 

Naturally, as the weight of the AZI contribution in- 
creases y^{B) behavior changes from decreasing to in- 
creasing. Note that the change of the sign of the de- 
pendence y^{B) has recently observed in a number of 
experiments. The proposed theory can, in principle, be 
useful in interpretation of some of these observations. 

In addition to this short discussion few comments are 
worth-wile, nevertheless, on the peculiar non-monotonic 
behavior of y{B) functions at small B, sometime with the 
small maximum at B = 0, which is found for some sets 
of parameters of models (see the curve 2 in Fig. 3c) . The 
fact is that such a behavior has recently been observed 
in some semiconductor devices fiS, In principle, the au- 
thors of this work assume that these non-monotonic de- 
pendences result from the exchange interaction between 
polarons at short distances. Our results show that some 
other mechanisms, for example the superposition of HFI 
and AZI induced relaxation, can lead to a similar y{B) 
behavior. 



4- Ag-mechanism 

The conventional Aij-mechanism implies the MFE gen- 
eration caused by quantum transitions between S and T 
states of the PP (|4.17p - (|4.19p . which result from the 
spitting Aw = {Ag)(3B of Zeeman frequencies (see Sec. 
IV.B)ji2 The difference of (7-factors of e- and /i-polarons 
results from the spin-orbital coupling of electron spins 
of polarons in semiconductors;^ i.e. the Ag'-mechanism 
can be considered as one of manifestations of the this 
coupling in MFEs. 

In our analysis of Ag-mechanism we neglect the effect 
of other (relaxation) mechanisms discussed above. 

The manifestation of Ag-mechanism in many types of 
MFEs is analyzed in detail in a large number of works 
as applied to different chemical and physical processes. 
The contribution of Ac^-mechanism to the PP recombina- 
tion yield ?/^(S) is also discussed in literature. Just this 
mechanism is believed to be responsible for y^(-B) \/B 



behavior found in some semiconductors at large fields 
B.-~!^ ' Moreover, in a number of works this behavior is 
used as a definition of the Ag-mechanism. Such a defini- 
tion is somewhat misleading. 

The fact is that Ag-mechanism, originally describing 
the effect of the term of the Hamiltonian H (|4.14p pro- 
portional to Ag, predicts different yf^{B) dependences, 
depending on the kinetics of relative motion of polarons: 

a. Free diffusion model. In the case of free diffu- 
sion the dependence yf{B) can approximately be rep- 
resented as^^i^ VfiB) ^ y/ Auj{l + ^jVAw)^^, where 
Auj = {Ag)[3B and ^ d/^/D. This formula shows for 
d^J A(jj/ D ^ 1 we get the dependence yf{B) ^ \J Auj ^ 
\fB (in some works considered as the manifestation Ag- 
mechanism), whereas in the opposite limit d\J Au / D > 1 
the dependence yf{B) saturates, i.e yf{B) ~ const. 

b. Diffusive cage model. In the diffusive cage model 
we consider the case of deep well, when the effect of recap- 
ture of escaping particles is small. In the deep well limit 
the term ^ fcg ~ •\/A in [eq. (|6.31l) ] can be neglected 
so that y^B) dependence (which becomes analytical) 
can approximated hy^^ y^{B) - {Auj f/[1 + ^^(Aw)^], 

where <^c ~ l|W^cll~^ l^^ weak reactivity limit S^c ~ 
W~^, as it is seen from eq. (|8.8p ]. This means that in the 
cage model Ag-mechanism predicts rapidly increasing be- 
havior of y^{B) at small B < \\Wj\/igl3): y^{B) - 
saturating at large B > |j Wg||/(g/3): y^B) ^ const. 

The simultaneous contributions of Ag-mechanism as 
well as HFI- and ASI-induced relaxation mechanisms can 
result in the additional specific features of the MFE shape 
y^{B). We are not going to discuss them in this work, 
but only mention that typically Ag-mechanism strongly 
contributes at large B, at which contributions of both 
other mechanisms are nearly independent of -B. In this 
case the change of MFE shape caused by Ag-mechanism 
can easily be identified and described if needed. 

IX. CONCLUDING REMARKS 

In this work PP recombination mechanisms of MFEs 
in disordered semiconductors are analyzed in detail. The 
magnetic field dependent PP recombination yield is dis- 
cussed as the most well known example of the MFE ob- 
servable. The hopping migration of polarons is assumed 
to result not only in the spatial evolution of polarons, 
but also in fluctuations of the HFI and AZI of polarons, 
thus leading to HFI and AZI induced spin relaxation. In 
our work we have considered the manifestation of these 
two mechanism of MFEs as well as Ag-mechanism. 

Simple analytical formulas are derived for the recom- 
bination yield Y{B) in two models of polaron migration: 
the free diffusion model and the model of diffusion in 
the well, or diffusive cage model. Analysis demonstrates 
that specific features of polaron migration shows itself in 
the shape of Y{B). Most important typical properties of 
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y(i?)-dependence for all three MFE mechanisms and in 
both models of migration are discussed in detail. 

Concluding this discussion we would like to point out 
two important points: 

1) Above-obtained results show that the width of 
the relaxation induced MFE (PP recombination yield) 
y^{B) is determined by the hopping rate: B^. ~ Wq/ {gl3), 

i.e. is much larger than typical HFI: B^ ^ y/ {B^). 

2) The proposed theory is applied to analysis of the 
magnetic field dependent yield of recombination oi e — h 
PPs, i.e. polarons with charges of opposite sign. It is 
clear, however, that obtained results are quite applicable 
to description of recombination of e— e PPs or h—h PPs 



as well. Just this kind of processes is recently considered 
as a mechanism of MFEs in organic semiconductors.'^^ 

3) This work concerns the analysis of recombination 
of polarons, i.e. particles with the electron spin 1/2. 
The obtained matrix formulas are, however, quite gen- 
eral and can be applied to describing similar processes 
with participation of particles with higher spins, for ex- 
ample, triplet exciton-polaron quenching,— li^i^i triplet- 
triplet annihilation, etc. The analysis of these pro- 
cesses is a subject of future works. 
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